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We revisit the Swift-Hohenberg model for two-dimensional hexagonal patterns in the bistability 
region where hexagons coexist with the uniform quiescent state. We both analyze the law of mo- 
tion of planar interfaces (separating hexagons and uniform regions), and the stability of localized 
structures. Interfaces exhibit properties analogous to that of interfaces in crystals, such as faceting, 
grooving and activated growth or "melting" . In the nonlinear regime, some spatially disordered het- 
erogeneous configurations do not evolve in time. Frozen states are essentially composed of extended 
polygonal domains of hexagons with pinned interfaces, that may coexist with isolated localized struc- 
tures randomly distributed in the quiescent background. Localized structures become metastable 
at the pinning/depiiming transition of interfaces. In some region of the parameter space, localized 
structures shrink meanwhile interfaces are still pinned. The region where localized structures have 
an infinite life-time is relatively limited. 

PACS numbers: 89.75.Kd, 05.45.-a, 64.60.My, 68.35.Dv 



I. INTRODUCTION 

Spatial disorder is a generic feature of pattern forma- 
tion. Whereas the simplest patterns that appear in num- 
bers of physical or chemical systems driven out of ther- 
modynamic equilibrium are periodic structures (such as 
stripes, squares or hexagons), in practice, many steady 
patterns observed in spatially extended systems are more 
complex and disordered at large scale Q, S I3- Why 
disordered configurations spontaneously form (without 
the help of impurities or any quenched disorder) and do 
not asymptotically evolve with time toward more ordered 
states is still not fully understood. 

For instance, systems where various (say, two) patterns 
of different symmetries coexist can easily lead to disor- 
dered states. A basic example is Rayleigh-Benard convec- 
tion of non-Boussinesq fluids, where hexagonal patterns 
bifurcate subcritically and are stable along with the uni- 
form conductive state within a whole interval of Rayleigh 
numbers. In bistability regions, asymptotic configura- 
tions may be either homogeneous (containing a single 
phase) or heterogeneous, with both phases present and 
distributed in a non-regular way. Heterogeneous states 
have been identified experimentally in Rayleigh-Benard 
convection 0,0, vibrated granular materials j^jOr Tur- 
ing patterns and related gas-discharge systems Spa- 
tial disorder has been also proposed as a generic feature 
of pattern formation in ecosystems, as illustrated by veg- 
etation patches in arid regions 0. 

On the theoretical side, the usual nonlinear model 
equations for periodic patterns often predict metastable 
disordered states as well. Regarding bistable/subcritical 
patterns, a lot of attention has been devoted in the past 
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years to localized structures, that are axisymmetric iso- 
lated structures of well-defined shape and size immersed 
in a uniform quie scent ba ckg round ( of vanishing order 
parameter) H M IH El El El [11 . Localized struc- 
tures have been identified as a possible origin of disorder 
in bistable systems, as asymptotic steady states are in 
some case composed of many localized structures ran- 
domly distributed in space d, . The stability of 
these "spatially chaotic" states, as opposed to periodic, 
can be explained by the particular nature of the effective 
pair interaction between two localized structures 0. 

In the present paper, we consider the Swift-Hohenberg 
model for hexagonal patterns in the bistability zone men- 
tioned above and study spatially more extended objects, 
namely, islands with hexagonal order surrounded by the 
quiescent background. A localized structure can be seen 
as a limiting case where such an domain is of very small 
size (of order Aq , the wavelength of the base pattern) . An 
interesting question to ask is up to what point localized 
structures and large islands differ from each other, al- 
though they probably share closely related features. We 
therefore analyze in the following the dynamics of a pla- 
nar interface separating an extended domain of hexagons 
and a uniform phase in the bistability region. Although 
interfaces have been relatively less studied than localized 
structures in this context |l3j |. we show that they play a 
central role in the large time evolution of two-dimensional 
extended systems, and in the formation of heterogeneous 
disordered states. 

After recalling a few known results on the weakly non- 
linear regime (Section^, we show that in the nonlinear 
regime interfaces can either move (with two modes: ex- 
panding or shrinking domains), or be pinned, depending 
on the value of the main control parameter (Section 
The theoretical treatment is based on an extension of 
the usual amplitude equation formalism to include "non- 
adiabatic" effects, according to some ideas presented in 
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El 111 111 E El. Some analogies of this system 
with interfaces in solid crystals are pointed out. The 
results are further confronted with numerical resolution 
of the Swift-Hohenberg equation in Section IIVI where a 
dynamical diagram for interfaces is drawn in parameter 
space. In Section we discuss the stability of hetero- 
geneous frozen states. We show that the stability region 
of localized structures only partially overlap the pinning 
regime of interfaces. We deduce that the existence of 
steady disordered configurations (and therefore of local- 
ized structures of infinite life-time) is possible in practice 
only in that pinning regime. Disordered frozen configu- 
rations obtained from random initial conditions can be of 
various types: they can contain localized structures, or 
localized structures and extended domains of hexagons 
("mixed" disorder), or extended domains and no local- 
ized structures. Changing the control parameters can 
drive the system outside of the pinning zone, keeping lo- 
calized structures stable: these are absorbed by growing 
domains and therefore become metastable. Conclusions 
are presented in Section IVTI 



II. MODEL AND KNOWN RESULTS ON 
INTERFACE MOTION 

We consider a system described by a dimensionless lo- 
cal order parameter ip(r,t), representing, for example, 
the mid-plane vertical velocity in convection problems. 
The spatio-temporal evolution of tp in non-Boussinesq 
convection can be described by a generalized Swift- 
Hohenberg equation 0| , 

^=6^-^(A + fcg)V + g2V'^-^^ (1) 

The parameter e ^ 1 is the main control parameter and 
represents the reduced Rayleigh number {R — Rc)/Rc, 
where R^ is the critical Rayleigh number, fcp = 27r/Ao 
is the characteristic wavenumber of the patterns, and 
the bare coherence length. In the case of convection, 
^0 is of order k^^ ^ its precise value depending on the 
boundary conditions at the fluid layer. The coefficient g2 
of the quadratic term in Eq.(^ models the strength of 
non-Boussinesq effects; 172 will be taken as positive in the 
following without restricting generality (172 <C 1). 

Equation or close variants of it, have been stud- 
ied in many contexts, for instance, to describe the 
microphase separation of asymmetric block-copolymer 
melts in the weak segregation limit |i23i | or for degener- 
ate optical parametric oscillators E3 ■ Recently, a similar 
equation was introduced to describe the evolution of veg- 
etation biomass (the local order parameter in that case) 
in arid climates [3]. 

Eq.JQ) can be written under a frequently used form: 

^ = Y.^-[^ + klf^ + G,^^-^\ (2) 
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FIG. 1: Schematic representation of an interface between 
a hexagonal pattern (left) and a uniform state (right). The 
colors roughly represent the sign of xp. The gray dots cor- 
respond to the next stable position of the interface in the 
pinning regime (see Section fill Bl belowl . 



with 

= (eo/2fco)^, T = ieo/Ot, 
S = e/(eo/4fco'), G2=g2/(eo/2fco). (3) 

We will take in the following e and ^0 (rather than e and 
(72) as the two main tunable parameters of Eq.(^|, 52 
being fixed to an arbitrary small value. The coherence 
length ^0 is thus considered as independent of Ao a priori. 
From (PI, decreasing ^0 is equivalent to increasing the 
bifurcation parameters E and G2. Keeping |e| ^ 1, the 
weakly nonlinear regime corresponds to ^0 ~ -^Oj while 
the nonlinear regime is reached with ^0 ^ -^o- 

For e < 0, the uniform solution -0 = (conductive 
state) is linearly stable, while it becomes unstable for 
e > 0. Standard weakly nonlinear analysis shows that a 
subcritical bifurcation to stable hexagonal patterns takes 
place at = — (^I/IS < 0. Hexagonal patterns remain 
stable up to e = 16g^/3, while stripes pattern become 
stable for e > 4(/|/3 [2J|. In the following, we restrict 
our study to the coexistence zone of the hexagonal and 
conductive states (em < e < 0). 

We come back in more details to a problem presented 
in fl3 | . We study a planar interface separating two semi- 
infinite domains, one composed of hexagonal patterns 
and the other in the uniform state ^ = 0. The local 
order parameter can be approximated at leading order 
(ei/2) as 

3 

V'o = ^ An{x, t) exp{iki ■ f) + c.c. , (4) 

n=l 

where the amplitudes An are slowly varying envelopes 
depending on cc, the coordinate along the axis x normal 
to the interface. The solutions of Eq.QJ satisfy the con- 
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dition 

lim An{x) = Ao > 0, lim An{x) =0 (n = 1, 2, 3). 

X — > — OO X — *oo 

The wave vectors fc„ are such that = ko, and make 
an angle of 27r/3 with each other. We choose ki with an 
orientation perpendicular to the interface: this situation 



corresponds a "faceted" hexagonal pattern having a raw 
parallel to the interface, as schematically displayed in 
FigH Faceted configurations are the more relevant ones, 
as further discussed in Section Hvl 

When |e| ^ 1, 32 ^ 1 and ^ ~ Aq, the equations 
satisfied by the amplitudes A„ can be derived from Eq. 
by standard multiple scale analysis 



dAi 
~W 
dA2 

dt 
dA3 

dt 



eA, + ^^dlA, + 252^2^ - H\Aif + 2\A2\^ + 2|^3|') 
eA2 + ^dlA2 + 252^1^ - 3(|^2p + 2\Ai\^ + 21^31') 
eAs + ^dlAs + 252^1^ - 3(|^3p + 2|Aip + 21^21'), 



(6) 
(7) 
(8) 
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{A is the complex conjugate of A). Note that A3 — A2 
from the planar geometry of the interface. 

We next recall some known results on the velocity v of 
the interface. Use the ansatz An = An{x — vt), multiply 
Eq^® (resp. Eqs.Q and ©) by dtAi (resp. dtA2 and 
df A^), add the equations and inte gra te their real part 
over X. The velocity is obtained as 



A3[45Ao/2-2g2] 
2D 



(9) 



where Aq is the amplitude of the regular hexagonal pat- 
tern 



and 



Ao^ 



D = 



92 + \j9l + 15e 
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(10) 



(11) 



The numerator of Eq. Q can be interpreted as a constant 
driving force acting on the interface, while D (which de- 
pends on the whole interface profile) plays the role of 
a friction term. The system of equations ©-((SJl has a 
potential form: there exists a functional "free-energy" 
F = J drT such that dtAn = -SF/SAn- The driving 
force in Eq.© (plotted below in Fig|21 as a function of 
e) is proportional to the difference of free-energy density 
T between the regular hexagonal state and the uniform 
state. Therefore, the sign of v is such that the domain 
with the lower free-energy expands at the expense of the 
other. From Eq.©, 172 being fixed, the velocity changes 
sign at e = eo = ~8g|/135 = (8/9)ein, a value located 
in the coexistence region. The uniform phase spreads if 
Em ^ e < eoi whereas the hexagonal phase does in the 
much wider interval < t <Q. The interface is station- 
ary only at e = eo. Therefore, except at the marginal 



value 7 the amplitude equation formalism at leading 
order predicts that spatially heterogeneous distributions 
of hexagonal and conductive regions should eventually 
evolve toward a single state, either hexagonal or conduc- 
tive, through interface motion. 



III. NON-ADIABATIC EFFECTS 
A. Modified amplitude equations 

From Eqs.©-®; the length scale of variation of the 
amplitudes, defining the interface thickness L, is roughly 
of order 



^o/v|e| or ^0/52- 



(12) 



Suppose now that |e| (or g|) is not that small compared to 
1, and/or ^0 is significantly smaller than Aq: the assump- 
tion i ^ Aq may not be fulfilled in that case. The spa- 
tial variations of An{x) are not necessarily much slower 
than the local oscillations of the patterns, and so-called 
"non-adiabatic" effects can appear. The effects of non- 
adiabaticity on moving interfaces have been studied the- 
oretically for one-dimensional patterns 0, l20 |. s tripe- 
stripe interfaces jU , square- hexagon interfaces |0| , and 
other cases . Non-adiabatic effects can be sufficient to 
pin interfaces and therefore to stabilize spatially hetero- 
geneous configurations of the order parameter. A weakly 
nonlinear analysis taking into account these effects in 
the present case is reported in Appendix ^ We de- 
rive the modified amplitude equations, with (small) non- 
adiabatic corrections that illustrate the coupling between 
the "fast" and "slow" length scales. Given the geometry 
depicted in FigQ] one obtains 
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ikQX 



+ id, [252(|Aip + + l^aP) - 6(Aa2A3 + A1A2A3)] 



fco 



9^2 

^3 = ^2- 



= -71- + (2ff2^2^i - A3 - SAjAs 



fco 



(13) 



(14) 
(15) 



In Eqs.(^-(^, the 5F/6Ai's are the right-hand sides 
of Eqs.©-©. 

The above equations involve directly the space variable 
X and therefore break the invariance by translation of 
standard amplitude equations. This is a consequence of 
the coupling between the amplitudes and the local phases 
of the waves. The additional terms in H13|l - H15(l can not 
be derived from a usual multiple scale analysis. Equa- 
tions lfT^ - ((TB|l have been derived in Appendix 1X1 in a 
lowest mode approximation (valid if L is still 3> Ao): 
the non-adiabatic terms oscillate in that case with the 
same wavenumber than the base pattern, fcp. Other 
terms oscillating as exp(2«fco), exp(3«fco)... should be kept 
if i ~ Ao, a situation corresponding to the strongly non- 
linear regime. 



B. Pinned interfaces for ^0 of the order of 

Fixing the non-Boussinesq parameter 52 to an arbi- 
trary small value, we determine in this Section the finite 



interval of values of the control parameter e such that 
interfaces do not move, or are pinned. This interval de- 
pends on the coherence length ^q, and can represent a 
significant part of the coexistence region when ^0 is small 
enough. 

To derive the modified law of motion of the interface, 
we use the method described in Section ^1 Denoting xq 
the position of the interface and using the ansatz An = 
An{x — xo{t)), one obtains from Eas. l(T^ - ((T3|l . 

w = — = j^ + j^ sm{koXo) (16) 



with 

F = Ai][45Ao/4-52], (18) 

and 



/oo 
dxcos(fcoa; -1- (p) [Aid, (852^? + 4^2^! - ISAiAj) 
-00 

+ 2^23,(252^1^2-3^3-3^2^2)], (19) 



D being given by (|ll(l . In Ea. ljlfill . the quantity p is a 
dimensionless positive number (p ^ 1) that can be iden- 
tified with the magnitude of a periodic pinning potential: 
the motion of the interface is analogous to that of a vis- 
cous particle sliding down over an inclined plane with un- 
dulations. From Eg. 1)16(1 . if p < the interface moves 
at speed v ^ 0, although non-constant in time. When 
p > \F\, the effective potential V seen by the interface 

has many local minima: there exist a set of discrete sta- 
(v) 

ble positions xj, for which v = 0. The interface is pinned 
and hexagonal domains do not grow nor shrink. Two suc- 
cessive steady positions are separated by a distance Ao. 
This value for the periodicity (result of a single mode 
approximation, see Appendix has a simple geometri- 
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cal interpretation: the closest steady configuration of the 
faceted hexagonal domain represented in Figure ^ is the 
one obtained by addition (or symmetrically, removal) of 
an interfacial "layer" of dots parallel to the interface (the 
gray dots of Fig^). This situation is analogous to that 
of a slow solidification (or melting) process controlled by 
activation barriers [26l |. 

In expressions (|11|) and (|19|) we have assumed that the 
AnS are real (their imaginary part is supposed to be 
small). To evaluate p, we solve numerically the first or- 
der equations (EJ-lCI) for the moving interface, and sub- 
stitute in H19|l . Although this solution is only approxi- 
mate (in the pinning regime, we expect some changes in 
the amplitude profiles due to the fact that interfaces are 
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FIG. 2: Strength of the driving force F and of tire pinning po- 
tential p (obtained for various values of the coherence length 
^o) as a function of t. (?2 = 0.3 here. 
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FIG. 3: Characteristic coherence length , defining the 
onset of a significant pinning regime, as a function of 32, as 
given by weakly nonlinear analysis. 



not moving), it captures the correct physical picture. A 
closer comparison with numerical solutions of Eq.JQ) is 
performed in Section ll VI below . 

In Figure 121 we plot the driving force \F\ and the pin- 
ning potential p given by Eas. lfTH|) - ((T^ as a function of e, 
for a fixed value of g2- We show some curves of p obtained 
for different values of the coherence length ^0 {F is in- 
dependent of ^o)- Decreasing ^0 increases non-adiabatic 
effects and the pinning potential: from Eq. H19|) , p roughly 
behaves as 



p ~ Aq exp(— aLfco), 



(20) 



with a a numerical constant of order unity and L given 
by lO. 

At fixed interfaces are pinned in a range einf(Co) ^ 
e < esup(^o)j while hexagons shrink for ei„ < e < einf 
and expand for esup < e, see Fig|2 For ^0 ^ (the 



theory 
numerics 



^„/X„(0.158 g,) 



FIG. 4: Relative magnitude of the pinning interval as a func- 
tion of a; = i^o/(0.158(?2Ao), as given by weakly nonlinear the- 
ory and numerical solution of Eq.Q. 



weakly nonlinear regime) , the size of the pinning interval 
is extremely small compared with |em|, and limited to the 
vicinity of eg ■ However, egup — Cinf blows up and covers a 
significant part of the coexistence region at a cross-over 
value i^\< Ao). The present calculation predicts that 
nearly the whole coexistence region (at least the interval 
[cq, 0]) has pinned interfaces for ^0 smaller than a critical 
value (see Fig|2l. We define the crossover value ^q^' as 
given by this critical value. In Figure |31 we plot ^q^'' as 
a function oi g2- A linear relation is observed: 



So 
Ao 



^ ~ 0.158 52. 



(21) 



Since 52 < 1, Co^'' < Ao in Eq.lEU- We conclude that 
pinning can be observed in practice only if the scale of 
variation of the amplitude is very short, corresponding 
to the nonlinear regime of Eq. , or equivalently if 



G2 > G 



(p) 



(22) 



in Eq.lPJ. Unfortunately, in that regime, the equations 
©-(HI used to derive these results are no longer valid: 
higher order terms should be included to improve accu- 
racy. We therefore do not expect this calculation to give 
a quantitative picture of pinning effects, but rather a 
qualitative one. An example of discrepancy between the- 
oretical and numerical results (presented in the following 
Section) is that the predicted pinning in the nearly whole 
coexistence region at low ^0 is not observed numerically. 
Instead, pinning only occurs in the lower fraction of the 
coexistence interval, as further shown in Sec. IIVI below. 

Nevertheless, the present weakly nonlinear analysis 
correctly captures the onset of a significant pinning 
regime for ^0 below a value , as well as the right 
order of magnitude for ^g^'' . In Figure 0] we plot the 
relative size of the pinning regime (csup — einf)/|em|, de- 
termined from theory and numerics, as a function of 
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X = ^o/(0-158(72Ao)- Due to the exponential form of 
Ea. H2(J|) . the theoretical curve is nearly a step function 
at X = 1. The numerical curve is softer, but with signifi- 
cant variations around a; = 1 as well. 



IV. NUMERICAL RESULTS 
A. Pinning of a planar interface 

We numerically solve Eq. Q by using a pseudo-spectral 
method and a time integration procedure described in 
[2^. The space is discretized on a square lattice, with 
2048 nodes along i, 128 nodes along y, and periodic 
boundary conditions. The initial condition is a rectan- 
gular hexagonal domain of length 1028 and width 128, 
comprised in-between two quiescent regions of -0 = 0. 
The two interfaces are faceted and parallel to the y di- 
rection. The lattice size Aa; set to unity. The base period 
Aq = 2n/kQ of the pattern is X. 

To determine whether interfaces are pinned or not, we 
compute the time evolution of the quantity 



dx 
Ao 



T ' 





(23) 



with and ly the box dimensions. a{t) is proportional 
to the length of the hexagonal domain, and therefore can 
be used to record the relative positions of the interfaces. 
Figure |31 shows transitions between moving and pinned 
boundaries as e is decreased, for a given 172 and ^o- As ex- 
pected, hexagonal domains are either expanding, steady 
or shrinking. In the moving regimes near a pinning tran- 
sition, the motion is noticeably non-uniform in time, in 
agreement with the form of the law of motion H16|l . 

Figure El displays a (e, ^o)-diagram representing the 
pinned and moving regimes of hexagon/uniform inter- 
faces. The parameter 32 is fixed to 0.2, but the shape 
of the diagram is independent of this value, which just 
sets the scales of both axis (gl for the x-axis and (72 for 
the y-axis, from Sections ITU and ITTi)) . For comparison, we 
also report the characteristic values of e given by weakly 
nonlinear analysis (^0 ~ Aq), introduced in Section ITTI 



and re-noted e 



(th) 



and 



here. When ^0 *C Ao, the 



coexistence region delimited by the lower value £„! be- 
comes wider as ^0 decreases: Cm deviates significantly 
from the value predicted by weakly nonlinear analysis. 
Yet, the picture drawn in the previous Section is qualita- 
tively reproduced. At any one can always identify two 
pinning-depinning transitions: For Em < e < einf, hexag- 
onal domains shrink; for einf < e < fsup, interfaces are 
pinned; for esup < e < 0, hexagons expand. For ^0 small 
enough, the pinning interval (csup — Cinf) quickly grows 
from almost zero and becomes comparable to |em|- The 
numerical cross-over value found numerically compares 
well with that given by formula H21|) , as shown in Figure 

The curve ein(Co) converges toward the value em^^ as 
^0 increases, while einf(^o) and esup(Co) converge toward 



1 r- 

0.98 - 

0.96 - 

0.94 - 

0.92 - 

a{t) 0.9 - 

o.ss - 

0.86 - 

0.84 - 

0.82 - 
0.8 



e = -0.0030 



e = -0.0035 

e = -0.0037 




(go = 0.2, fo = 0.03Ao) 



n 500 1000 1.500 2000 2500 3000 3500 4000 4500 5000 



t 



FIG. 5; Position of an interface as a function of time, ob- 
tained by numerical solution of Eq.Q, for various values 
of the control parameter e. The two upper curves display 
a net motion over pinning barriers, the hexagonal domain 
expanding into the uniform region (e > esup — —0.0037). 
The two following curves correspond to pinned interfaces 
(tint < e < Esup), and the lower one to shrinking hexagonal 



domains (em < e < einf). 



Jth) 



On the other hand, e^n and e-mf tend to merge as 
Co/Ao — > 0, a trend also reproduced by the theoretical 
curves of Fig (21 

In the region located below the dotted line of Figure El 
(low ^0 and |e|), fronts are moving, with possible instabil- 
ities toward labyrinthine patterns or non-trivial uniform 
phases. This region is not of interest here and was not 
investigated. 



B. Faceting of oblique interfaces 



One can show from Section HIH that interfaces have to 
be faceted to be pinned. As detailed in the Appendix, 
if none of the wavenumbers ki of the base pattern is di- 
rected along the normal to the interface, the nonadiabatic 
terms (IA2p will always integrate to zero along the trans- 
verse coordinate y. In p^ . it was similarly argued that 
oblique interfaces could not be pinned. We show here 
that oblique interfaces actually modify their shapes to 
become pinned. More generally, in the pinning regime, 
an hexagonal domain surrounded by a conductive state 
"0 = evolves toward a polygonal shape compatible with 
its three wavenumber directions. This situation can be 
seen as a nonequilibrium analogue of faceting in crystals 
due to strong anisotropy in surface tension |27l| . 

As an example, we have solved numerically Eq. JQ) with 
an initial condition composed of a planar interface with 
a wavenumber fci making an angle 9 = 5.73° with the in- 
terface normal (Figure [7^). The interface evolves toward 



7 




0.01 - 



-0.008 -0.007 -0.006 -0.00.5 -0.004 -0.003 -0.002 -0.001 

e 

FIG. 6: Dynamics of hexagons/uniform interfaces in param- 
eter space (e, i^o), for g2 = 0.2. (The shape of the diagram 
does not depend on the value of g2, that only sets the scales 
of the axis.) The limiting values eJJ and eo'' given by weakly 
nonlinear analysis (^o/^o ~ 1) are reported for comparison. 
The region below the dotted curve is not of interest here, see 
text for details. 




FIG. 7; Faceting of an oblique interface in the pinning regime 
(p2 = 0.2, e = -0.0041, ^0 = O.OSAo). The initial interface 
makes an angle 8 — 5.73° with the normal, a.) t — 100; b) 
t = 1900. 



a wavy, faceted shape (Figure 03). 



V. DISORDERED STATES: COMPETITION 
BETWEEN DOMAINS OF HEXAGONAL ORDER 
AND LOCALIZED STRUCTURES. 

From the preceding analysis, we conclude that spatially 
disordered configurations of the order parameter can be 
stabilized via pinning effects as soon as the coherence 
length is of order ^q^'' (or G2 in Eq. Q of order , given 
by H22|l l. These configurations can contain immobile, ex- 
tended hexagonal domains with faceted boundaries. In 
this subsection, we wish to compare the previously deter- 



mined dynamical regimes for planar interfaces with the 
stability diagram of localized structures. 

Localized structures are axisymmetric (dot-like) struc- 
tures of well defined shape and size frequently encoun- 
tered in bistability regions. They have been widely stud- 
ied and have been recognized to play an irnportant role 
in the production of disordered states 0, Hini. The 
stability of assemblies of localized structures irregularly 
distributed in space is relatively well understood. The ef- 
fective interaction force between two localized structures 
changes sign and oscillates with their separation distance 
(a property reminiscent of Ea. (|16|l for interfaces), there- 
fore allowing many stable discrete positions that are un- 
likely to lead to periodic arrangements of structures . 
We show below that the likelihood of observing local- 
ized structures in disordered solutions of Eq.l^ in fact 
largely depends on whether planar boundaries are pinned 
or not. We find that localized structures can actually co- 
exist with pinned extended domains. In other cases, they 
are metastable, i.e. susceptible to be swept by any ex- 
panding domain of hexagonal order. 

We have determined the stability diagram for a single 
localized structure. Calculations were performed on a 
square of 128 x 128 grid points over an extended period of 
time (20000 time units), and with a small additive noise 
term in Eq.l^. The initial localized structure (a region 
of size ~ Ao/2 with -0 > 0) can either shrink, converge 
to a stable form, or expand. When localized structures 
expand, they generate hexagonal order, as described in 
details in refs. 0,0|- Similarly to extended domains, 
we find that localized structures shrink if e < e[^f^\ are 

stable if e[^f^ < e < esup\ and expand if esup^ < e. 
The results are reported in the (^qj e)-plane in Figure |H1 
together with the data of FigureElfor planar fronts. Note 
that the stability zone of localized structures is located 
toward the right of the pinning regime of interfaces. 

The stability interval eiup — ^inf^' and the pinning in- 
terval Esup — Cinf for boundaries become significant in the 
same range of values of the coherence length ^o- The 
hatched region of Fig|Sl indicates where both intervals 
overlap. In this region of parameters, localized struc- 
tures actually have an infinite life-time. Anywhere out- 
side this zone, localized structure are either unstable, or 
metastable, i.e. susceptible to be swept by a moving front 
of expanding hexagons. Also notice that when ^0 /-^o is of 
order 1, eiup\— ^inf^) tends to 0~. Therefore, locahzed 
structures always shrink in the weakly nonlinear regime. 

To further illustrate these results, we have solved the 
Swift-Hohenberg equation |^ with 172 = 0.2 and random 
initial conditions, in various regions of Figure|Sl For con- 
venience, we fix ^o(= O.O3IA0) in the vicinity of \ The 
initial condition for is given by a random Gaussian field 
of zero mean and variance A^, with a a dimcnsionless 
tunable parameter of order 10. 

Inside the hatched region of Figure|Hl(e|jff'^'' < e < Csup) 
and with e relatively closer to ^\f^f\ the asymptotic pat- 
terns obtained are mainly composed of aperiodic assem- 
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FIG. 8: Stability diagram of localized structures (LS), superposed to that of Figure El for planar fronts {g2 — 0.2). In the 
hatched region, fronts are pinned and localized structures do not expand nor shrink, and therefore have an infinite life-time. 




FIG. 9: Different kinds of frozen configurations of -i/; (in gray 
scale) obtained at large times from random initial conditions 
((72 = 0.2). a) The coherence length is = O.OSlAo and 

( LS) 

e — —0.0055 is near tj^j in the hatched region of Figure 
IHl (t = 20000, a = 13). b) Same ^o, with e = -0.0040, 
closer to esup in the hatched region of Figure |S|(t — 240000, 
a = 8.66). The pinned hexagonal clusters are faceted, and 
stable localized structures are still observed, c) = 0.04Ao 



and 



Einf 



< 



-0.0045 < Min[e „ , Esup] : interfaces are 
pinned but localized structures are unstable (shrink). The 
smallest stable structures are clusters with two dots. 



biles of localized structures of liquid-like character, with 
no visible hexagonal order. The example shown in Figure 
IHK at e = 0.0055 (for a system of size 256 x 256) is similar 
to the patterns presented in various studies on "spatial 
chaos" 0, 0, 131 . Given a random initial condition, at 
short times, most of the pattern quickly relax to the con- 
ductive state tp = 0, with only a few remaining patches 
that further lead to immobile localized structures. One 
can control the final density of localized structures by 
means of the parameter a in the initial condition. The 
configurations become more dilute as a decreases. 



If e is increased and set closer to esup m \eLt '^e^, 



ALS) 

^ini ' "=supj, 

random initial conditions gradually tends to form other 
kinds of frozen disordered states. These states are 
"mixed" , characterized by the presence of both local- 
ized structures and faceted extended hexagonal domains 
(Figure IHti). The interfaces in Fig. are pinned, 

as expected. The enhancement of hexagonal order for 
quenches closer to egup can be qualitatively understood 
from Section Fill Bl the driving force F increases with e 
faster than the pinning potential p, even in the pinning 
zone, sec Fig|5] The time evolution of a typical configu- 
ration is the following: A few stable localized structures 
form at intermediate times. Small clusters, composed 
of 2,3,4 or more dots, are also present. Clusters with 2 
dots do not grow (nor move). Some clusters with 3 or 
more dots are unstable: some of them, that are initially 
irregular, expand slowly by the creation of successive lay- 
ers of dots (they may absorb in the process stable local- 
ized structures standing in their neighborhood) . Growing 
clusters eventually reach steady states, corresponding to 
faceted polygons with an underlying hexagonal structure. 
These extended hexagonal domains (generalizations of 



FIG. 10: Pattern evolution for esup < e = —0.0026 < eLp • (The other parameters are ~ O.OSlAo and (?2 = 0.2.) Four 
successive configurations of a same run are shown, respectively at times t = 5000, 10000, 15000, 20000. Localized structures 
are stable, but eventually absorbed by the growing domains. 




FIG. 11: Same parameters as in Figure ITUl with e = —0.0010 > eiup^ ■ Three successive configurations at times t — 
2000, 3000, 4000 are displayed. Transient localized structures are now unstable, and nucleate fast-growing hexagonal domains. 



the 7-dots structures observed in ref.[l3|) coexist with 
the remaining single localized structures. It takes a long 
time (typically 300, 000 times units) for the whole system 
to reach a steady state. The normal to an interface can 
be directed along non-crystalline directions, due to, for 
instance, merging processes taking place between neigh- 
boring clusters. Domains can be fairly extended by the 
time they reach a faceted shape, that further prevent 
their growth. 

The facets ensure stability, although some defects 
(penta-hepta dislocations) can be present in the bulk of a 
domain. A slight cusp inward is visible on an interface in 
the upper left corner of FigEb- This feature is reminis- 
cent of a grooving effect with sharp angles, a phenomenon 
well known in polycrystals with a strong anisotropy in 
surface tension |28l |. 

Further increase in e beyond Cgup leads, as expected, 
to homogeneous states of hexagons, with no conductive 
regions at large times. Figures [TUl shows successive con- 
figurations obtained at intermediate times in the case 
Esup < e(= —0.0026) < eiup^ Small hexagonal clus- 
ters grow, in a way analogous to a stable phase in a 
metastable one. Although the single localized structures 
present are individually stable, they are absorbed dur- 
ing growth, as hexagon/uniform interfaces sweep them 



in a finite time. Hence, localized structures can not be 
observed at arbitrary large times and are metastable in 
this range. The "absolute" stability region of localized 
structures is in practice restricted to the hatched zone of 
FiglSJ included in the pinning zone of interfaces. 

(LS) 

When the control parameter exceeds esup , transient 
short-lived localized structures may still emerge from ran- 
dom initial conditions. They further destabilize by form- 
ing structures of six-fold symmetry and then growing 
bubbles, as illustrated in the sequence of Figure [TTl Lo- 
calized structures nucleate hexagonal order in this case, 
as already observed in the Lengyel- Epstein model of Tur- 
ing patterns |13| and experimentally in semiconductor- 
gas discharge systems 0- 

An other possible case, illustrated by Figure Et, cor- 
responds to pinned interfaces and unstable (shrinking) 
localized structures. It happens in the roughly triangu- 

(LS) 

lar region delimited by the curves einf, Csup and e^^f 
in Fig|S| In FigO:, the smallest structures observed are 
clusters with two dots. Finally, for e < Cinf, any initial 
random condition converges toward the conductive state 
ij = 0. 
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VI. CONCLUSIONS 

We have studied the dynamics of interfaces in the 
bistabihty region of conductive and hexagonal states of 
the Swift-Hohenberg equation We have examined 
the consequences of the various dynamical regimes of 
interfaces on the occurrence of spatial complexity and 
on the metastability of smaller structures like localized 
structures. Steady disordered states are possible in the 
nonlinear regime only, and for a choice of parameters such 
that planar interfaces remain pinned. In the expansion 
regime (e > Csup) localized structure may be stable but 
are susceptible to be swept by any traveling interface in 
a finite time. 

The size of the pinning interval [emfjCsup] is signifi- 
cant as soon as the bare coherence length is of the 
order of (or smaller than) a characteristic value ^q^"* , ap- 
proximately given by formula H21(l . An extended domain 
with pinned interfaces can be seen as a nonequilibrium 
analogue of a solid crystal. Domains can contain dis- 
locations in the bulk, and even vacancies. Defects can 
induce grooving-like effects at the interface (non-planar 
interfaces with a cusp inward). 

On general grounds, decreasing the coherence length 
(keeping the other parameters constant) increases nonlin- 
ear effects and enhances hexagonal order (V'(t = oo) 7^ 0), 
as illustrated by the monotonic curves of em (Co) and 
einf(Co) in the diagram of FigEl The quantity can 
be seen as a measure of the pattern "strength" , a "weak" 
pattern yielding the phase ip = 0. Though, one of the 
nontrivial features of Figl^is the non-monotonic behav- 
ior of Csup with Co- There exist a small interval of e's 
such that pinned states undergo a depinning transition 
(the system becoming completely filled with hexagons) if 
one slightly increases Co- In this case, a "weaker" pattern 
fills space better. 

Heterogeneous frozen configurations can be of various 
kinds: they can include localized structures ( "liquid" - 
like disorder), or both localized structures and extended 
faceted domains of the condensed phase ( "mixed" disor- 
der); or even, faceted domains with small clusters and 
no isolated localized structures. As the Swift-Hohenberg 
equation ^ is the minimal model allowing subcritical 
patterns and bistability, we expect the results presented 
here, specially the conclusions drawn from the diagrams 
of Figsini and |S| to be quite general and applicable to 
a variety of systems. Let us mention Turing chemi- 
cal patterns, nonlinear optical systems and related de- 
vices 0, 01- (Rayleigh-Benard convection and block- 
copolymer melts are unfortunately less promising candi- 
dates as their Co is relatively large |1,|23|.) Our results 
show in particular that crystal-like patterns can be con- 

I 



^ „ dx' fy+^y dy\ , 

Jx ^0 Jy ly 



tained inside polygonal domains compatible with sym- 
metries, and of arbitrary shape. Such patterns could 
be used as lithographic templates for technological ap- 
plications (like micro-electronics) that require patterned 
surfaces over well-defined regions of space. 
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APPENDIX A: AMPLITUDE EQUATIONS WITH 
NON-ADIABATIC EFFECTS 

The equations satisfied by the amplitudes An can be 
derived from Eq.(^ by standard multiple scale analy- 
sis The spatial variations of ip^ see Eq.Q), involve 
two kinds of length scales: a short scale associated with 
the periodicity of modulations, Ao, and a large one as- 
sociated with the scale of variation of the amplitudes, 
of order L (the interface thickness). When e ^ and 
52 — >■ 0, L ^ Ao, and it is possible to decompose spatial 
variations into "fast" variables, f = (x,?/), and a "slow" 
variable X = e^^^x (as well a a slow temporal variable 
T = et), such that A„ = A„{X,T). The equations for 
the A'^s are derived from a solubility condition for at 
the higher order following i/jq. The solubility condition 
for the amplitude An reads, 

rX + \0 J / fly J I ^ 

/ ^ lim / ^[i„(^o)+52V'o'-V'o']e-''="'^ = 0, 

Jx '^O 'b^°oJo h 

(Al) 

where the linear operator L„ = e — ec^y — eCo/(4fco)^[2(A;n- 
x)dxdxY follows from Eq.lQ), V'o being given by Q and 
kn — kn/\kn\. In (|A1|I . we have taken the limit ly — *■ 
00 for commodity, since the interface profile is invariant 
along y. Close to onset, X and x can be considered 
as two independent variables and the only non- vanishing 
contributions of the integral (|A1I) come from the terms 
oscillating as exp(ifc„ • r ') in the brackets. Coming back 
to physical variables, this leads to the well-known coupled 
amplitude equations (|SJ)-||HJ) 

According to relation \i'2.l . when Co is smaller than Aq 
(or when e and 172 are not small compared to unity), the 
"fast" and "slow" variables [x and X, respectively) can 
not be considered as independent in the solubility condi- 
tion (|A1|I for An- ipo being given by the decomposition 
if one expands, for instance, the nonlinear term 32 "00 
m some following integrals appear: 



')A„,{x ') exp{iK^x ') eyipiiKyy ') (A2) 
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with the notation A-i = Ai and —3 < m, Z < 3. In 
Ea. HA2|l . K = K^x + Kyij = ki + fc,„ — fc„, with the 
notation k-i — —ki. This time, let us consider the non- 
resonant terms only, i.e. those with K ^ 0. The An 
profiles of the planar interface depending on x ' only, the 
contribution / is exactly zero unless Ky = 0. If Ky = 
and Kx ^ 0, then / reduces to 

dx' 

I^2g2 —Ai{x')Am{x')exp{iKxx'), (A3) 

Jx ^0 

where is actually a multiple (positive or negative) 
of fco- The usual amplitude formalism assumes that Ai 
varies slowly, so that Ai{x') =cst in the interval [x,x + 
Ao]. In that case, the non- resonant / given by Ea. ljA3(l 
always vanishes. However, as soon as L/Xq is large but 
finite, the complex exponential does not exactly integrate 
to zero if one sits in the interface region; /, although 
small, does not vanish. 

Anticipating on Section llllBl fsee Ea l20() . one can show 
that a term like (|A3() in the solubility condition (|A1() gives 
an additional contribution of order 

^exp{~a\Kx\L)smiKxXo) (A4) 

to the velocity v of the boundary given by Eq.©. (See 
further Eq.(Iin|)-(IIll).) In Ea. (|X4)l . a is a constant of or- 
der unity, xo the position of the interface with respect to 
an arbitrary origin, and D given by (|11|) . \Kx\ can take 
the values fco, 2ko, 3ko... In HA4|I . \Kx\L is finite but still 
assumed to be ^ 1: one can keep among the nonadia- 



batic terms / those with the smallest \Kx\ {\Kx\ = ko), 
and neglect those involving higher wavenumbers. 

One now derives the amplitude equations with their 
non-adiabatic corrections. One develops 52V'o ~ V'o (with 
ipo given Eq.Q), and, for each solubility condition HA1|) 
(n — 1,2,3), one has to find out, with the notation 
introduced above, which are the products AiAm (and 
AjAiAjn) such that K = ki+km — kn (or kj+ki+km — kn) 
has Ky = and = ±A:o. Note that one must elimi- 
nate among the eligible products those that are resonant 
in any of the leading order amplitude equations 10-® or 
their complex conjugates. These resonant terms, added 
to similar ones coming from the linear operator L„ in 
(jAl|) . constitute one of the 6 standard amplitude equa- 
tions: the sum of these terms cancels out at order e'^/^ 
(the order of the terms in the brackets of Eg. HAlj) '). and 
thus would contribute in non-adiabatic effects at the fol- 
lowing order, not considered here. For example, in de- 
riving the nonadiabatic contributions of the equation for 
A2, one finds an eligible term (|A2I) involving a product 
A1A2, characterized by \K\ = \ki\ = feg. But this term 
is resonant in the equation for A3, and hence needs to 
be discarded. The products of amplitudes appearing in 
the leading nonadiabatic terms can not be already known 
products. 

In order to calculate the integrals (|A3p . we next assume 
that the amplitudes vary slowly in the interval [x, x + 
Ao] and An{x') ~ An{x) + {x' — x)dxAn{x). After some 
tedious but straightforward algebra, equations l|13|l - l|15|) 
are obtained. 
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